Unconventional superconducting pairing symmetry induced by phonons 
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The possibility of non-s-wave superconductivity induced by phonons is investigated using a simple 
model that is inspired by Sr2Ru04. The model assumes a two-dimensional electronic structure, a 
two-dimensional spin-fluctuation spectrum, and three-dimensional electron-phonon coupling. Taken 
separately, each interaction favors formation of spin-singlet pairs (of s symmetry for the phonon in- 
teraction and d x 2_ y 2 symmetry for the spin interaction), but in combination, a variety of more 
unusual singlet and triplet states are found, depending on the interaction parameters. This may 
have important implications for Sr2Ru04, providing a plausible explanation of how the observed 
spin fluctuations, which clearly favor d x 2_ y 2 pairing, may still be instrumental in creating a super- 
conducting state with a different (e.g., p-wave) symmetry. It also suggests an interpretation of the 
large isotope effect observed in Sr2Ru04. These results indicate that phonons could play a key role 
in establishing the order-parameter symmetry in S^RuCU, and possibly in other unconventional 
superconductors. 

PACS numbers: 74.20.Rp, 74.20.Mn, 74.70.Pq 



I. INTRODUCTION 

The existence of superconducting states that break 
symmetries of the normal state beyond gauge symme- 
try has generated considerable interest. Such uncon- 
ventional pairing was first discovered in 3 He, in which 
the three superfluid phases all have spin-triplet (S = 1) 
atomic Cooper pairs with relative orbital angular mo- 
mentum L = 1 (p-wave)i Strong evidence for uncon- 
ventional superconductivity was subsequently found in 
heavy-fermion materials^ In the high-T c cuprates, phase- 
sensitive experiments have provided convincing evidence 
for L = 2 (d-wave) spin-singlet pairing^ The pairing 
symmetry in Sr2Ru04, a non-cuprate perovskite oxide 
superconductor, has been unambiguously shown to be 
lower than the symmetry of the underlying crystal lat- 
tice, i.e., L > 0. 4 It is generally believed that Sr 2 RuC>4 
is a triplet p-wave superconductor? 4 , though it has been 
noted that a <i-wave singlet state is also possible^ It has 
been suggested that some compounds could be supercon- 
ductors with L > More exotic types of pairing, such 
as order parameters with odd-frequency dependence and 
mixed singlet-triplet Cooper pairs, have also been con- 
sidered for real materials iSiS 

By itself, the electron-phonon interaction is always at- 
tractive, and hence it maximizes the pairing energy for 
an order parameter that is positive everywhere on the 
Fermi surface. Conventional wisdom thus suggests that 
L > pairing symmetries must be induced by non- 
phonon mechanisms, such as magnetic or Coulomb in- 
teractions. A number of authors have pointed out, how- 
ever, that while phonons are always more pairing in the 
s channel, they may also be pairing in other channels as 
well^fi If, for example, strong on-site Coulomb repulsion 
suppresses s-wave pairing because of the large overlap of 
the wave functions of the paired particles, or if there is an 
additional pairing contribution in an L > channel from 



other (non-phonon) mechanisms, unconventional pairing 
that results mainly from electron-phonon coupling is pos- 
sible. Less well appreciated is the idea that a combination 
of electron-phonon coupling and another pairing interac- 
tion (e.g., spin-fluctuation exchange) may result in a su- 
perconducting state that is not the ground state of either 
of the two interactions taken separately. 

In this paper we consider a simple model, loosely based 
on Sr2Ru04, that illustrates this effect. The model as- 
sumes a highly two-dimensional electronic structure, and 
it includes a two-dimensional spin-fluctuation-induced 
interaction that favors d-wave pairing of x 2 — y 2 (or 
cos 2</>) symmetryjAi and a three-dimensional electron- 
phonon interaction that favors an s-wave state. The 
spin interaction is based on the experimentally observed 
spin- fluctuation spectrum of Sr2Ru04, which has negli- 
gible dispersion along While less is known about 
the electron-phonon coupling in Sr2Ru04, the three- 
dimensional nature of the crystal structure suggests that 
the electron-phonon interaction in the material may be 
much more three dimensional than the electronic struc- 
ture itself. When the coupling strengths of the two in- 
teractions are comparable, we find that the symmetry 
of the emerging ground state can be changed by tun- 
ing the widths of the interactions. Among the states 
that are found to be stable in certain regions of phase 
space are two-dimensional (z-independent) p-wave states 
and three-dimensional (z-dependent) singlet and triplet 
states, all of which arise from the interplay between the 
two interactions. 



II. MODEL 

The spin dependence of the pairing in a superconduc- 
tor can be taken into account using a 2 x 2 matrix A(k) 
to represent the gap function. Within this spin matrix 
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formalism, the BCS gap equation has the form 
A SlS2 (k)= £ K 2 ™(k,kOA S3S4 (k')*\k',T), (1) 

where 

V SlS2S3Si (k,k') = (-ksi;-ks 2 \V\ -k's^k'ss) (2) 

is the matrix element of the effective electron-electron 
pairing interaction V, for which we adopt a convention 
in which positive V corresponds to an attractive inter- 
action, and F is a function of the temperature T and 
the quasiparticle energy E^m^ In a crystal with inver- 
sion symmetry, the pairing potential can be separated 
into spin-triplet (+) and spin-singlet (-) channels, 

i/ slS2S3S4 (k,k') = y(+)(k,k')-sw s3S4 +y(-)(k,k')5(-) 2S3 

(3) 

where 

'-'.S1S2S3S4 = (^S1S4^S2S3 ^ ^S 1 S 3 S S2 s 4 ) /2. (4) 

The gap equation decouples into separate equations for 
the singlet and triplet channels. For singlet pairing, the 
gap matrix has the form &o(k)i<j y , while for triplet pair- 
ing, it has the form [a ■ d(k)]ia y , where a are Pauli spin 
matrices. Since the total wave function must be anti- 
symmetric under interchange of particles, the scalar or- 
der parameter Ao(k) for singlet pairing must be even in 
k and the vector order parameter d(k) must have odd 
parity. 

Near T = T c , the gap function is small and the gap 
equation can be linearized and reformulated as an eigen- 
value problem: 

/(T c )A slS2 (k)= ]T ^ 2SlS3S4 (k,k')A S3S4 (k'). (5) 

k'S3S4 

In the singlet channel, this becomes an eigenvalue prob- 
lem for y(~'(k, k'), with eigenfunctions given by the 
scalar order parameter Ao(k). For triplet pairing, the 
eigenvectors of V^ +S) (k, k') correspond to components of 
the vector order parameter d(k). Among eigenfunctions 
that have spatial symmetry compatible with the spin 
pairing, the one with the largest eigenvalue f(T c ) pro- 
duces the first superconducting instability that manifests 
as temperature is lowered, and hence defines the physical 
T c . 

In our model, the pairing potential has two contribu- 
tions: an electron-phonon interaction V p h that is always 
attractive regardless of the symmetry of the spin pair- 
ing, and an interaction V sp that is mediated by anti- 
ferromagnetic spin fluctuations. The latter is attractive 
in the triplet spin channel but repulsive in the singlet 
spin channel. Therefore, the total pairing interaction is 
V = (V ph + V V )SM + (V ph - V sp )S(~\ 

The model system considered here is inspired by 
Sr2Ru04. A tetragonal, three-dimensional crystal struc- 
ture with Z?4/j symmetry is assumed, with lattice con- 
stants a and c. The electronic structure is highly two- 
dimensional, with a cylindrical Fermi surface of radius 



k F = 0.97r/a, similar to the 7 sheet of the Fermi sur- 
face in Sr2Ru04. The electron-phonon interaction is 
assumed to be three dimensional in nature while the 
antifcrromagnetic-spin-fluctuation-induced interaction is 
two dimensional. Both interactions depend only on 
q = k — k'. The spin-fluctuation-induced interaction 
V sp (q) is taken to be an isotropic Gaussian in q x and 
q y , centered on q = (27r/3a, 27r/3a, q z ) and equiva- 
lent lines in the Brillouin zone [i.e., (27r/3a, 47r/3a, q z ), 
(4-7r/3a, 27r/3a, q z ), (47r/3a, 47r/3a, q z )], in accordance 
with both theorjii^ and experiment^ The electron- 
phonon interaction V p h{<\) is a three-dimensional Gaus- 
sian function centered at (n/a, 7r/a, 7r/c), i.e., the corners 
of the Brillouin zone. The electron-phonon interaction is 
allowed to be anisotropic, with a different width in the 
, z direction than in the plane. The anisotropy of the 
electron-phonon interaction is characterized by the pa- 
rameter a — aw ph /cWp h , where the widths are defined in 
terms of the Gaussian form exp[— (q/w) 2 ]. 

Because the Fermi surface in the model is perfectly 
cylindrical and the interaction potentials are separable in 
Cartesian coordinates, the eigenvalue problem separates 
into z and planar parts, with product eigenfunctions. 
The in-plane eigenfunctions are of the form cos(m(/>) 
and sin(rn0), while the z eigenfunctions are of the form 
cos(nfc 2 c) and sin(nfc z c), where m and n are integers and 
tan</> = k y /k x . For n ^ 0, the two-dimensional spin 
interaction cancels out since V sp (k, k') cos(nk' z c) — 
Sk' Vs P (k, k') sin(nfc^c) = if V sp has no z dependence. 
This means there are degenerate spin-triplet and spin- 
singlet states that have the same in-plane eigenfunction. 
The symmetry of the in-plane eigenfunction determines 
which one has the cosine solution and which one has the 
sine solution in k z . In a real material, this degeneracy 
between triplet and singlet states would be broken by a 
variety of effects, including warping of the Fermi surface 
and the spin-orbit interaction. 

Numerical calculations were carried out using a dis- 
crete grid of k points on the Fermi surface. The circum- 
ference of the Fermi surface was sampled with 40 equally 
spaced grid points, while 20 points were used along k z . 
The grid was selected so that the spacing between grid 
points was small compared to the the widths of the inter- 
actions. We have also checked convergence by increasing 
the number of grid points to 60 x 30 and 80 x 40. 

The effective strength of each interaction in coupling 
electrons on the Fermi surface is measured by 

^sp(ph) = ^ V sp ( ph )(k,k'). (6) 
k,k'GFS 

For fixed ratios of R = X p h/^sp, phase diagrams were 
mapped out as a function of interaction widths. Because 
of the discrete sampling of k-points, results are presented 
only for interactions with Gaussian width greater that 
0.l7r/a in the plane. To construct the phase diagrams, 
we used a recursive algorithm that takes advantage of 
the fact that most phases occupy large regions in phase 
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space. The gap equation is solved for all points along the 
border of the rectangular phase diagram. If the ground 
state is the same at all points on the border, all points in 
the interior are assumed to have the same ground state. 
Otherwise, the region is divided into two rectangles and 
the process is repeated recursively. For each phase di- 
agram, phase space was sampled with a range of reso- 
lutions, from 5 x 5 to 400 x 400 points, to mitigate the 
chance of missing phases that might be stable in only 
small interior regions of the diagram. 



III. RESULTS 

We consider first the case when one of the interactions 
dominates. When the spin coupling is much larger than 
the phonon coupling (R <C 1), only one phase appears 
regardless of the widths of the spin and phonon interac- 
tions. This is the singlet cos 2(f) state (commonly denoted 
k% — ky). Because the spin interaction is repulsive in the 
singlet channel, to induce pair formation, it must scatter 
electrons between points k and k' on the Fermi surface 
that have order parameter of opposite sign (see Eq. EJ. 
For cos 2(f> symmetry, three of the four lines in the Bril- 
louin zone on which the spin interaction is maximal do 
this; the fourth couples nodal regions. Hence this d-wave 
state is stabilized by the antiferromagnetic spin fluctua- 
tions. When the phonon interaction dominates (i? > 5), 
only an s-wave spin singlet appears in the phase diagram. 
This is expected since the attractive electron-phonon in- 
teraction is pair forming when the order parameter has 
the same sign at points k and k' on the Fermi surface. 

Figure ^ shows the phase diagram for equal spin and 
phonon coupling strengths (R = 1). In the diagram, 
the width of the spin-fluctuation-induced interaction in- 
creases from left to right, and the width of the electron- 
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FIG. 1: Phase diagram for coupling-strength ratio R — 
Aph/A S p = 1 and electron-phonon anisotropy a — 1. See 
Table Q] for the symmetry of states that appear in this phase 
diagram. 



TABLE I: Symmetry of gap functions of states that appear in 
Figs. 1-5. For orbital symmetry (L), the commonly used s-, p- 
, d-wave, etc. nomenclature is adopted, with s* , for example, 
denoting extended s-wave symmetry. For spin-singlet states 
(S — 0), basis functions for the scalar order parameter Ao are 
given, while for spin triplets (S = 1), basis functions for the 
components of the vector order parameter di (i = x, y, z) are 
listed. In the weak-coupling, non-relativistic approximation 
used in this work, all triplet states with the same k depen- 
dence are degenerate, independent of the direction of the d 
vector. For instance, all 11 p-wave states listed for tetrago- 
nal crystals in the Sigrist and Ueda review" (zk x , itk x ±yk y , 
etc.) are degenerate in this approximation and correspond 
to the state labeled here as A. Actual gap functions include 
higher-order Fourier components that transform in the same 
way under the symmetry operations of the crystal (D^) as 
the basis functions shown in the table. 

S L Ap(k) or dj(k) 

Alp sin (f), cos (f> 

B s 1 

C d cos 2(f) 

D d* cos 2</>cos(fc z c) 

1 / cos 2(f)sm(k z c) 

E s* cos(2fc z c) 

1 p* sin(2fc z c) 

F d sin 2(f) 

G g sin 40 

H d sin (f>sm(k z c), cos (f> sin(fc z c) 

1 p* sin </>cos(fc z c), cos<?icos(fc z c) 

I d* cos 2(f) cos(3fc z c) 

1 /* cos 2(f) sin(3fc z c) 

J s* cos(4fc z c) 

1 p* sin(4fe z c) 

K d* sin 0sin(3fc z c), cos (f> sin(3fc z c) 

1 p* sin 0cos(3fc z c), cos (f> cos(3fc z c) 

L d* sin 2d> cos(2fc z c) 



phonon interaction increases from top to bottom. The 
anisotropy parameter for the width of the electron- 
phonon interaction is a = 1. In the regions labeled 
B (green), the s-wave singlet favored by the electron- 
phonon interaction is stable, while in region C (blue), the 
cos 2(p spin-singlet state favored by the spin-fluctuation- 
induced interaction is stable. New phases also appear. 
There is a manifold of degenerate triplet states with in- 
plane p-wave pairing (A, orange), and degenerate singlet 
and triplet states that have ^-dependent order parame- 
ters (D, purple). The character of these phases is summa- 
rized in TableQ] We have mapped out phase diagrams for 
coupling-strength ratios ranging from R = to R = 5. 
With the electron-phonon anisotropy fixed at a = 1, no 
states other than the ones that are shown in the R = 1 
diagram are found to be stable. 

It is interesting that the two-dimensional triplet man- 
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ifold (A, orange) emerges despite the fact that the in- 
teractions individually stabilize singlet states. For this 
triplet state, the coupling matrix (Eq. EJ) has two degen- 
erate eigenvalues, with eigenfunctions cos <f> and sin <fi, for 
each component of the vector order parameter d. This 
generates the set of allowed p-wave states on a cylindri- 
cal Fermi surface with tetr agon al lattice symmetry, as 
derived, for example, in Ref. I 1. Our model corresponds 
to the weak-coupling, non-relativistic regime, where these 
states are all degenerate. A chiral member of this set of 
p-wave triplet states, with d oc {k x + ik y )z, is a lead- 
ing contender for the order parameter in Si^RuO^ The 
two-dimensional triplet manifold is stable over a range 
of coupling strength ratios, but in all cases, it is limited 
to the region of the phase diagram corresponding to nar- 
row spin interactions. This can be understood from the 
form of the spin interaction. In the triplet channel, the 
spin-fluctuation-mediated interaction is attractive, so to 
induce pairing it must scatter electrons between points on 
the Fermi surface that have order parameters of the same 
sign. Wave vectors along the line q = (47r/3a, 47r/3a, q z ), 
for which the spin interaction is strong, have an in-plane 
component that is slightly larger than 2kp, the diame- 
ter of the Fermi cylinder. Hence, if the spin interaction 
is narrowly peaked in the planar directions, these wave 
vectors contribute little to the total spin coupling. With 
wider spin interactions, however, they have a strong pair- 
breaking effect if the order parameter has p-wave symme- 
try in the plane. 

The singlet phase in region D (purple) has symme- 
try cos2(/>cos(fc z c) and can be thought of as an extended 
d-wave state, as illustrated in Fig. [21 It belongs to the 
same representation as the purely two-dimensional cos 2<p 
state. For small w p h, it is more stable than the cos2(/> 
state. This is because in the two-dimensional cos 2<f> state 
favored by the spin interaction, the phonon interaction 
centered at q = ("ft/a, ir/a, ir/c) is always pair breaking, 
while the inclusion of a cos(fc z c) factor allows the phonon 
interaction to become pair forming. Degenerate with this 
extended d-wave singlet is a set of triplet /-states, one 
example of which is d oc cos20sin(fc z c)z. Note that be- 
cause the two-dimensional spin interaction integrates out, 
the stabilization of the extended d-wave singlet and the 
degenerate /-wave triplets is due entirely to the electron- 
phonon interaction. 

Figure [3] takes vertical and horizontal cuts through 
the phase diagram shown in Fig. ^ to illustrate how 
the eigenvalues change with the widths of the phonon 
and spin interactions. In Fig. [3^, the width of the 
electron-phonon interaction is fixed at w^ h = 0.6257r/a 
with a = 1. In Fig. |2t>, the width of the spin- fluctuation- 
induced interaction is fixed at w sp — 0.257r/a. In addition 
to the four phases that appear in Fig. corresponding 
to those states with the largest eigenvalues within the 
phase space of the diagram, phases with other symme- 
tries appear in Fig. [3] Table [I] lists the symmetry of 
these states. Curve H (gold), for example, corresponds 
to a manifold of three-dimensional singlets with Aq given 




FIG. 2: Order parameter of the singlet state stable in region D 
(purple). White corresponds to the maximum positive value 
of Ao and black corresponds to the minimum negative value. 
The F point is at the center of the cylinder. The full eigen- 
function for R = 1, w sp = 0.57r/a, w^ h — 0.57r/a, and a — 1 
is plotted, but the lowest-order non-zero Fourier component, 
cos2<^cos(fc z c), clearly dominates. 



by linear combinations of sin 0sin(fc z c) and cos 0sin(fc z c), 
which includes the chiral d-wave state, (k x + ik y ) sin(/c z c) 
that is compatible with existing experimental data on the 
symmetry of the order parameter in S^RuO^ When 
both the spin and phonon interactions are very narrow, 
or when both are wide, this three-dimensional d-wave sin- 
glet (H) becomes more stable than the two-dimensional p- 
wave triplet (A) that is generally believed to be the most 
likely candidate for the pairing symmetry in Sr2Ru04. 4 
The increased stability of the d-wave singlet as compared 
to the p-wave triplet in certain parameter regimes arises 
from the same mechanism that was discussed earlier in 
relation to the extended d-wave state (D, purple). Other 
symmetries, however, can make even better use of the 
structure of the pairing interaction in these parts of pa- 
rameter space, so the three-dimensional d-wave singlet 
does not become the favored phase. 

By narrowing the width of the electron-phonon inter- 
action in the z direction, the stabilizing mechanism for 
states with fc z -dependent order parameters can be en- 
hanced. Figure 01 illustrates how the eigenvalues change 
with anisotropy in the phonon interaction. The eigen- 
values depend on a only if the order parameter varies 
with k z . As expected, the three-dimensional d-wave sin- 
glet (H, gold) is further stabilized as a is increased, but 
within our model and within the parameter space ex- 
plored, it never emerges as the most stable phase. Nev- 
ertheless, one can imagine that with a different structure 
of electron-phonon coupling in reciprocal space, this so- 
lution could become stable in some parameter regime. 
When Wp h is very small, an extended s-wave singlet state 
(E, dark red) with symmetry cos(2fc z c) becomes stable 
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FIG. 3: Dependence of eigenvalues on (a) width of spin- 
fluctuation- induced interaction, w 3p , and (b) in- plane width 
of electron- phonon interaction, w^ h , for coupling-strength ra- 
tio R — 1 and electron-phonon anisotropy a = 1. In (a), 
w lh = 0.6257r/a, and in (b), w sp = 0.25ir/a. Dashed lines 
indicate states with fc 2 -dependent eigenfunctions. See Table 
[I] for the full symmetry of the eigenfunctions. 



FIG. 4: Dependence of eigenvalues on the anisotropy of 
the electron-phonon interaction, for coupling-strength ratio 
R — 1, spin-interaction width w 3p = 0.57r/a, and in-plane 
phonon-interaction width w ph = 0.6257r/a. Dashed lines in- 
dicate states with fc z -dependent eigenfunctions. See Table 
for the full symmetry of the eigenfunctions. 
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FIG. 5: Phase diagram for coupling-strength ratio R — 1 
and electron-phonon anisotropy a = 2.5. See Table for the 
symmetry of states that appear in this phase diagram. 



IV. CONCLUSIONS 



and occupies a large area of the phase diagram (Fig. [SJ. 
As discussed earlier, for states with fc z -dependent order 
parameters, the two-dimensional spin fluctuations inte- 
grate out and play no role (i.e., they are neither pairing 
nor pair breaking) . This is reflected in Fig. [5] by the 
fact the the phase boundary between region D (purple) 
and region E (dark red) is perfectly horizontal. 16 Each 
of those phases represents one singlet (extended d or ex- 
tended s, for D and E, respectively) and several triplet (/ 
or extended p for D and E, respectively) states. Similar 
/ states have been discussed in connection with heavy- 
fermion superconductors, e.g. in UPt3^£ 



The combination of spin- and phonon-mediated pairing 
interactions with different structure in momentum space 
can lead to a rich variety of unconventional pairing states, 
many with symmetries and parities that could never be 
stabilized by either interaction alone. Within the toy 
model considered here, with a cylindrical Fermi surface, 
two-dimensional spin-fluctuation-mediated pairing, and 
three-dimensional phonon-mediated pairing, there are 
two-dimensional pairing states with no dispersion in z 
and three-dimensional fc 2 -dependent states. For the two- 
dimensional states, the interplay of the phonon and spin 
interactions, each of which favors a two-dimensional sin- 
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glet state, can result in a triplet state. For the three- 
dimensional states, the symmetry of the model sup- 
presses the effect of the spin interaction, and phonon- 
induced non-s-wave states become possible. 

This toy model is only loosely connected with the un- 
conventional superconductor S^RuCU, yet it suggests a 
resolution to the seeming paradox regarding the sym- 
metry of the superconducting state in this compound. 
While most researchers believe that superconductivity 
is induced by the strong spin fluctuations observed in 
this system, the momentum dependence of these fluc- 
tuations, well-known from neutron measurements 
should stabilize the same d x 2_ y 2 pairing symmetry as in 
the cuprates. Yet this symmetry can be excluded based 
on a variety of experiments^^ On the other hand, a 
large isotope effect on T c indicates the importance of 
phonons^ yet the s or extended s symmetry favored by 
phonons can be excluded with even more confidence. Our 
results suggest the unconventional and rather unusual (in 



all likelihood, chiral) symmetry of the superconducting 
state in Sr 2 Ru04 emerges as a result of interplay of elec- 
tron coupling with both phonons and spin fluctuations, 
although neither of these agents separately can stabilize 
such a state. A more quantitative application of these 
ideas to Sr2Ru04 requires a complete mapping of the 
electron-phonon coupling in the Brillouin zone to com- 
plement the existing map of spin fluctuations, as well as 
consideration of the full, multi-sheeted Fermi surface. 
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